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We discuss the status of the solutions of the solar neutrino problem in terms of oscillations of v e into active 
or sterile neutrinos. We present the results of a global fit to the full data set corresponding to the latest Super- 
Kamiokande (SK) data on the total event rate, their day-night dependence and the recoil electron energy spectrum, 
together with the data from Chlorine and Gallium experiments presented at the J/-2000 conference. We show the 
possible solutions in the full parameter space for oscillations including both MSW and vacuum, as well as quasi- 
vacuum oscillations (QVO) and matter effects for mixing angles in the second octant (the so called dark side). 
We quantify our results in terms of allowed regions as well as the goodness of the fit (GOF) for the different 
allowed solutions. Our conclusion is that from the statistical point of view, all solutions for oscillations into active 
neutrinos: the large mixing angle (LMA), the low mass (LOW), the small mixing angle (SMA) and the QVO 
solutions are acceptable since they all provide a reasonable GOF to the full data set. The same holds for the 
SMA solution for oscillations into sterile neutrinos. LMA and LOW-QVO solutions for oscillations into active 
neutrinos seem slightly favoured over SMA solutions for oscillations into active or sterile neutrinos but these last 
two are not ruled out. We also analyze the dependence of these conclusions on the uncertainty of the SSM 8 B flux 
and on the removal of the data from one of the experimental rates. We also present the results in the framework 
of four neutrino oscillations which allows for oscillations into a state which is a combination of active and sterile 
neutrino. 
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1. Introduction: Experimental Status 

The sun is a source of v' e s which are produced in 
the different nuclear reactions taking place in its 
interior. Along this talk we will use the v e fluxes 
from Bahcall-Pinsonneault calculations [1] which 
we refer to as the solar standard model (SSM). 

At present these neutrinos have been detected 
at the Earth by seven experiments which use dif- 
ferent detection techniques: 
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Homestake (chlorine) [2] 

SAGE [3] and GALLEX+GNO [4,5] (gallium) 



71 



Ga -> 71 Ge + e" 



- Kamiokande and SK [6,7] v e e scattering on wa- 
ter 

Due to the different energy threshold for the de- 
tection reactions, these experiments are sensitive 
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to different parts of the solar neutrino spectrum 
as represented in Fig. 1. They all observe a deficit 
between 30 and 60 % which seems to be energy 
dependent mainly due to the lower Chlorine rate. 
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Figure 3. Recoil electron energy spectrum mea- 
sured by SK normalized to the SSM prediction. 



Figure 2. Measured solar neutrino event rates 
normalized to the SSM prediction. 



To the measurements of these six experiments 
we have to add also the new results from SNO [8] 
first presented in this conference. They are how- 
ever still not in the form of definite measured 
rates which could be included in this analysis. 

In this conference SK has also presented their 
results after 1117 days of data taking on: 

- The recoil electron energy spectrum: SK has 
measured the dependence of the even rates on 
the recoil electron energy spectrum divided in 
18 bins starting at 5.5 MeV. They have also re- 
ported the results of a lower energy bin 5 MeV 
< E e <5.5 MeV, but its systematic errors are 
still under study and it is not included in their 
nor our analysis. The spectrum shows no clear 
distortion with xjiat = 13/(17efof). 17 dof=18 
bins-1 free normalization. 

-The zenith Angle Distribution (Day/Night Ef- 
fect) which measures the effect of the Earth Mat- 
ter in the neutrino propagation. We have in- 
cluded in the analysis the experimental results 
from SK on the zenith angle distribution of events 
taken on 5 night periods and the day averaged 
value. SK finds few more events at night than 
during the day but the corresponding Day-Night 
asymmetry 
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Figure 4. Zenith angle distribution measured by 
SK normalized to the SSM prediction. 



In order to combine both the Day-Night in- 
formation and the spectral data SK has also 
presented separately the measured recoil energy 
spectrum during the day and during the night. 
This will be referred in the following as the day- 
night spectra data which contains 2 x 18 data 
bins. 

The most generic and popular explanation of 
the solar neutrino anomaly is in terms of neutrino 
masses and mixing leading to oscillations of v e 
into an active [y^ and/or i/ T ) or sterile neutrino, 
v s . We pass now to discuss some issues that have 
been raised lately in the literature concerning the 
computation of the corresponding neutrino sur- 



vival probabilities in the full range of mass and 
mixing relevant to the solar neutrino problem. 

2. Survival Probabilities: QVO and the 
Dark side 

The presence of neutrino mass and mixing im- 
ply the possibility of neutrino oscillations. For 
solar neutrinos we know that depending on the 
range of mass and mixing v e can undergo oscil- 
lations cither in vacuum [9] or via the matter- 
enhanced MSW mechanism [10]. However this 
broken picture of solar neutrino oscillations con- 
tains a set of approximations which, as we discuss 
next, are not longer phenomcnologically valid 
when performing the analysis of the solar neu- 
trino data. In order to clarify this issue let us 
first review the calculation of the solar neutrino 
survival probability in the two-neutrino mixing 
case. 

The survival amplitude for a solar v e neutrino 
of energy £ at a detector in the Earth can be 
written as: 

2 

A ee - ]T A s el Af e exp[-im 2 (L _ r)/2E] (2) 

Here A^, L is the amplitude of the transition u e — ► 
Vi (yi is the i-mass eigenstate) from the produc- 
tion point to the Sun surface, Af e is the ampli- 
tude of the transition Vi — ► v e from the Earth 
surface to the detector, and the propagation in 
vacuum from the Sun to the surface of the Earth 
is given by the exponential. L is the distance 
between the center of the Sun and the surface 
of the Earth, and r is the distance between the 
neutrino production point and the surface of the 
Sun. The corresponding survival probability P ee 
is then given by: 



P ee = Pi Pie + P 2 P 2e + 2v/PiP 2 Pi e P 2e cos£ (3) 

Here Pi = |^4fJ 2 is the probability that the so- 
lar neutrinos reach the surface of the Sun as |fj), 
while Pi e = |^4fe| 2 is the probability of Vi arriving 
at the surface of the Earth to be detected asac e . 
Unitarity implies Pi + P 2 = 1 and P le + P 2e = 1. 
The phase £ is given by 
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where S contains the phases due to propagation 
in the Sun and in the Earth and can be safely ne- 
glected. In the evaluation of both Pi and P 2e 
the effect of coherent forward interaction with 
the Sun and Earth matter must be taken into 
account. 

From Eq. (3) one can recover more familiar ex- 
pressions for P ee : 

(1) For Am 2 /E < 5 x KT 17 eV, the matter 
effect suppresses flavour transitions both in the 
Sun and the Earth. Consequently, the probabili- 
ties Pi and P 2e are simply the projections of the 
v e state onto the mass eigenstates: Pi = cos 2 9, 
P 2e = sin 2 9. In this case we are left with the 
standard vacuum oscillation formula: 



py 



1 - sin 2 29 sin 2 (Am 2 (L - r)/4E) 



(5) 



which describes the oscillations on the way from 
the surface of the Sun to the surface of the Earth. 
The probability is symmetric under the change 
of octant 9 «-► -| — 9 and change of mass sign 



Am 2 



■Am . Notice that the simultaneous 



application of both symmetries translates simply 
into a relabelling of the mass eigenstates v\ <-» v 2 . 
Therefore only one of these symmetries is physi- 
cally independent. This means that we can take 
Am 2 > without loss of generality and keep 
< 9 < 77. Then the symmetry of the survival 
probability under the change of octant 9 <-> f — # 
implies that each point in the (Am 2 , sin (29)) pa- 
rameter space corresponds to two physically inde- 
pendent solutions one in each octant. 

Averaging Eq.(5) over the Earth Orbit L(t) = 
Lq[1 — ecos27r^] one gets : 
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where e = 0.0167 is the orbit eccentricity and Jo 
is the Bessel function. In Fig. 5. a we display the 
value of (P e y e ac ) as a function of AE/Am 2 . As seen 
in the figure for large values of Am 2 the probabil- 
ity averages out to a constant value 1 — | sin 2 (29). 

(2) For Am 2 /E > 10~ 14 eV, the last term in 
Eq. (3) vanishes and we recover the incoherent 
MSW survival probability. In this case Pi and 
P 2e must be obtained by solving the evolution 
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Figure 5. P ee as a function of 4E/Ar 



equation of the neutrino states in the Sun and 
the Earth matter respectively: 
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e + ^cos2<9 -^sin20 
^sin 2(9 V x - ^ir cos 20 




(7) 



with 

v » = v ^-\^ 



K, = o 



(8) 



where N e i n \ (r) is the electron (nucleon) num- 
ber density which are proportional to the Sun or 
Earth matter density. In Fig. 6 we show N e (r) 
for the Sun. As seen in the figure the Sun density 
profile for r < O.9i?0 can be very well approxi- 
mated by an exponential N(r) = Nq exp(— r/rO) 
with constant exponent slope r$ = i? Q /10.54. 
The approximate solution for the evolution in 
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Figure 6. Solar density profile for BP2000 model. 



the Sun takes the well-known form 



Pi = 2 + (5 - P e) cos(2^ m>0 ) 



(9) 



where P c denotes the standard Landau-Zener 
probability [11] and # m .o is the mixing angle in 
matter at the neutrino production point: 
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with Ao = 2E(V e — Vx) evaluated at the produc- 
tion point, c2 = cos(26*) and s2 — sin(2(9). For 
the approximation of exponential density profile 

A 2 

7 = 7r " r which is independent of the point 
in the Sun where the resonance takes place. Im- 
provement over this "constant slope" exponential 
density profile approximation can be obtained by 
numerically deriving the exact N e (r) profile at the 
resonant point. In this case 7 = 7r-^-r (r res ). 

The physical interpretation behind Eq.(9) is 
very simple. At the production point v e has a 
projection over the mass eigenstate v\ given by 
the mixing angle in matter 6 m fi. This neutrino 
evolves adiabatically till the resonant point where 



both mass eigenvalues become closer and the neu- 
trino has a probability P c of "jumping" into the 
other mass eigenstate (or staying in the state V\ 
with a probability 1-P C ). During the day no Earth 
matter effect is to be included and the survival 
probabilities P- ie are obtained by simple projec- 
tion of the v e state onto the mass eigenstates 
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2e.DAY 



= 1 — Pie. day = sin 9 and one obtains 



Pie.DAY = o + (o - P c) ^s(29 mfi ) COs(2tf) (12) 



In Fig. 5.b ( 5.c) we plot the this survival proba- 
bility as a function of 4E/Am 2 for large (small) 
mixing angle. 

Let us make some remarks concerning Eq.(12): 

(i) In both limits of very large and very small 
E/Am 2 pMSW _> j _ i sin 2(26») (s ee Fig.5.b). 
For large E/Am 2 , P c = cos 2 9 and cos{29 m fi) = 
-1. While for small E/Am 2 , P c = and 
cos(26' m .o) = cos(20). So in principle one could 
expect a perfect connection between the asymp- 
totic small E/m 2 probability for the vacuum os- 
cillation regime and the large E/m 2 behaviour of 
the MSW regime. However as seen when com- 
paring panels (a) and (b) of Fig. 5 there is an in- 
termediate range, 2 x 10 14 < AE/Am 2 < 10 16 
eV, where both adiabaticity is violated and the 
cos£ coherent term should be taken into account. 
The result is similar to vacuum oscillations but 
with small matter corrections. We define this 
case as QVO [14-16]. The range of E/m 2 for 
the QVO regime depends on the value of E/m 2 
for which the MSW probability in Eq. 12 acquires 
the asymptotic value 1 — | sin 2 (29), the smaller 
E/m 2 the more separated the MSW and vacuum 
regimes are, and the narrower the QVO region is. 

(ii) Due to matter effects P<^ sw is only 
symmetric under simultaneous (Am 2 , 9) — > 
(—Am 2 , | — 9). For Am 2 > the resonance is 
only possible for 9 < j and MSW solutions are 
usually plotted in the (Am 2 , sin 2 (20)) plane as- 
suming that now each point on this plane rep- 
resents only one physical solution with 9 in the 
first octant. But in principle non- resonant solu- 
tions are also possible for 9 > j, the so called 
dark side [18-20,22]. 

It is clear from these considerations that in or- 
der to compute the survival probability for so- 



lar neutrinos, valid for any value of the neutrino 
mass and mixing, the full expression (3) has to 
be evaluated. The results presented in the follow- 
ing sections have been obtained using the general 
expression for the survival probability in Eq. (3) 
with Pi and P2e found by numerically solving the 
evolution equation in the Sun and the Earth mat- 
ter. For Pi we use the electron number density of 
BP2000 model [12]. For Pi e we integrate numer- 
ically the evolution equation in the Earth matter 
using the Earth density profile given in the Pre- 
liminary Reference Earth Model (PREM) [13]. 

Before moving to the next section we want to 
discuss the validity of the analytical and semi- 
analytical approximations based on the exponen- 
tial profile for the sun density in the evaluation 
of the survival probability when used to describe 
the matter effects in the QVO regime as well as 
for mixing angles in the second octant. In order 
to illustrate this point we display in Fig. 7 the 
survival probability P ee (without the oscillating 
cos ^ term) as a function of 4E/Am 2 for different 
values of the mixing angle obtained by our nu- 
merical solution as well as from the corresponding 
analytical approximations (12). We discuss two 
analytical approximations: the exponential ap- 
proximation with constant slope tq — _Rq/1054 
and with point dependent slope ro(r). 

We see from Fig. 7 that the analytical results 
differ from the results of the numerical calcula- 
tions. In particular, the analytical result with 
constant exponential slope shows a "slower" tran- 
sition to the vacuum oscillation regime or, in 
other words, it overestimates the size of the mat- 
ter effects in the QVO region. The same value of 
the survival probability appears for about twice 
larger E/Am 2 . This is due to the fact that for 
such values of E/Am 2 the adiabaticity breaking 
occurs very close to the Sun edge where the den- 
sity falls much faster than the exponential ap- 
proximation as shown in Fig. 6. Similar con- 
clusions have been drawn in Refs. [15-17]. On 
the other hand the analytical approximation with 
point dependent slope shows a "faster" transi- 
tion to the vacuum oscillation regime or, in other 
words, it underestimates the size of the matter 
effects in the QVO region. 

Notice also that, originally, Eq. (11) was de- 
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Figure 7. Survival probabilities as a function of 
4E/Am 2 for different values of the mixing angle. 
The full line corresponds to the numerical solu- 
tion while the dashed (dotted) lines corresponds 
to the analytical approximation for exponential 
density profile with constant (point dependent) 
slope. 



the "pure" MSW regime 4E/Am 2 < 10" 14 eV" 1 , 
while they are very small for QVO. 

(b) Earth Regeneration is always positive (it 
always leads to an increase of the v flux at night) 
for mixing angles both in the first or in the second 
octant. 
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Figure 8. Regeneration of v & in the Earth matter 
as a function of the AE/Am for different mixing 
angles. 



rived for a mixing angle 6 < ? where resonant 
enhancement is possible. However for the con- 
stant slope exponential profile, 7 is constant (in- 
dependent of the resonant condition) and both 
Eq. (10) and Eq. (11) can be analytically contin- 
ued into the second octant, 9 > 7r/4, and used 
to compute the corresponding survival probabil- 
ity. This is illustrated in Fig. 7. As seen in the 
figure, for values of the mixing angle close to max- 
imal mixing — 7r/4, the survival probability is 
mirror-symmetric in the first and second octant, 
while, as expected, the symmetry breaks down 
as we depart from maximal mixing. We also see 
that the analytical approximation overestimates 
the size of matter effects in the QVO region on 
the second octant as well. 

Finally we want to comment briefly on the 
Earth matter effects. In Fig. 8 we show the Earth 
regeneration factor Re = P2c,night - P2c,day 
(obtained numerically) as a function of AE/Am 2 
for two values of the mixing angle in the first and 
second octant. We see that: 

(a) Earth Matter effects are only relevant in 



3. Two neutrino Oscillations 

We now describe the results of the analysis of 
the solar neutrino data in terms of v e oscillations 
into active or sterile neutrinos. For details on the 
statistical analysis applied to the different observ- 
ables we refer to Ref. [21]. 

We first determine the allowed range of oscilla- 
tion parameters using only the total event rates of 
the Chlorine, Gallium and SK experiments shown 
in Fig. 2. For the Gallium experiments we have 
used the weighted average of the results from 
GALLEX+GNO and SAGE detectors and we do 
not include the Kamiokande data as it is well in 
agreement with the results from the SK experi- 
ment and the precision of this last one is much 
higher. 

Using the predicted fluxes from the BP98 
model the \ 2 for the total event rates is Xssm = 
60 for 3 d. o. f. This means that the SSM together 
with the SM of particle interactions can explain 



the observed data with a probability lower than 

io- 12 . 

The allowed regions in the oscillation param- 
eter space are shown in Fig. 9. In the case of 
active-active neutrino oscillations we find that 
the best-fit point is obtained for the SMA so- 
lution which has a probability or goodness of the 
fit of 50 the LMA and LOW solutions. In the re- 
gion of vacuum oscillations we find a "tower" of 
regions around local minima corresponding to the 
oscillation wavelength being an entire fraction of 
the Sun-Earth distance. In Table 1 we give the 
values of the parameters in these minima as well 
as the GOF corresponding to each solution. 

Notice that following the standard procedure, 
the allowed regions for a given set of observables 
are defined in terms of shifts of the x 2 function for 
those observables with respect to the global min- 
imum in the plane. Defined this way the size of 
a region depends on the relative quality of its lo- 
cal minimum with respect to the global minimum 
but from the size of the region we cannot infer the 
actual absolute quality of the description in each 
region. In order to give this information we list in 
Table 1 the GOF for each solution obtained from 
the value of x 2 at the different minima. We see 
that for oscillations into active neutrinos the best 
solution by far for the description of the rates is 
the SMA.' 

For oscillations into sterile neutrinos we find 
only the SMA solution and some small region 
for LOW-QVO. The LMA and most of the LOW 
solutions are not acceptable for oscillation into 
sterile neutrinos. Unlike active neutrinos which 
lead to events in the SK detector by interacting 
via neutral current (NC) with the electrons, ster- 
ile neutrinos do not contribute to the SK event 
rates. Therefore a larger survival probability for 
S B neutrinos is needed to accommodate the mea- 
sured rate. As a consequence a larger contribu- 
tion from S B neutrinos to the Chlorine and Gal- 
lium experiments is expected, so that the small 
measured rate in Chlorine can only be accommo- 
dated if no Be neutrinos are present in the flux. 
This is only possible in the SMA solution region, 
since in the LMA and LOW regions the suppres- 
sion of 7 Be neutrinos is not enough. Notice also 
the SMA region for oscillations into sterile neu- 



trinos is slightly shifted downwards as compared 
with the active case. This is due to the small 
modification in the neutrino survival probability 
induced by the different matter potentials. The 
matter potential for sterile neutrinos is smaller 
than for active neutrinos due to the negative NC 
contribution proportional to the neutron abun- 
dance (see Eq. (8)). For this reason the resonant 
condition for sterile neutrinos is achieved at lower 
Am 2 . 

On the other hand SK data on the day-night 
variation of the event rates and the spectrum 
data lead to excluded regions in the parameter 
space. Since both the observed zenith angular 
dependence and the recoiled energy spectrum are 
compatible with the no-oscillation hypothesis (up 
to an overall normalization which is taken free), 
these measurements translate in the exclusion of 
those regions of the oscillation parameter space 
where the expectation for those observables is far 
different from the no-oscillation one. In this way 
we see in Figs. 10, 11 and 12 the regions excluded 
by the zenith angle, average spectrum and com- 
bined day-night spectra data. The zenith angular 
data excludes the region where large regeneration 
in the Earth is expected (See Fig. 8) which for 
neutrino energies of few MeV (above SK thresh- 
old) occurs for few 10~ 8 < Am 2 /eV 2 < 10~ 5 . 
In Table 1 we show the effect of the inclusion of 
this observation in the position of the different lo- 
cal minima. We see that when combined with the 
zenith angle data the LMA minimum is shifted up 
and the LOW minimum is shifted down towards 
values of parameters where the Earth regenera- 
tion is smaller. 

The excluded region from the non observation 
of any clear distortion of the energy spectrum is 
shown in Fig. 11. Its main effect is to suppress 
the vacuum oscillation solution where a large dis- 
tortion of the spectrum is expected due to the 
"imprints" of the oscillation wavelength. On the 
other hand for the SMA a positive tilt of the spec- 
trum is also predicted due to the raising with the 
neutrino energy of the survival probability (see 
Fig. 5.c). For this reason the local SMA mini- 
mum shifts to lower values of the mixing angle 
where a flatter probability is expected. On the 
other hand both LMA and LOW solutions pre- 
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Figure 9. 90, 95 and 99 % CL allowed regions from the analysis of the total rates only. The global 
minimum is marked with a star while the local minima arc denoted with a dot. 
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Figure 10. 99 (shadow), 95 and 90 excluded regions from the analysis of SK measured zenith angle 
dependence. The inverted triangle represents the minimum used to define the regions. 
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Figure 11. 99 (shadow), 95 and 90 excluded regions from the analysis of the dally averaged SK measured 
spectrum. The inverted triangle represents the minimum used to define the regions. 



diet a rather flat spectrum and in consequence 
the position of these minima are little affected. 

Finally Fig. 12 shows the excluded region from 
the SK day-night spectra which contains simul- 
taneously the information on the Earth regen- 
eration effect and the energy dependence of the 
survival probability. Combining statistically this 
information with the data from the total rates 
(what implies 3+2x18 data points with 1 free nor- 
malization for the spectra) we obtain the allowed 
regions which we display in Fig. 13 while the cor- 
responding GOF for the global solutions are given 
in Table 1. There are some points concerning 
these results that we would like to stress: 

(a) Despite giving a worse fit to the observed 
total rates, once the day-night spectra data is in- 
cluded the LMA gives the best fit. This is mainly 
driven by the flatness of the spectrum and it was 
already the case with the last year data as pointed 
out inRef. [21]. 

(b) The GOF of the LOW solution has in- 
creased considerably as it describes the spectrum 
data very well despite it gives a bad fit to the 
global rates. LOW and QVO regions are con- 



nected at the 99 %CL and they extend into the 
second octant so maximal mixing is allowed at 99 
% CL for Am 2 in the LOW-QVO region. Notice 
that part of this allowed region would be missing 
when making the old fashioned analysis in terms 
of pure MSW or vacuum oscillations. 

(c) What is the situation for the SMA solu- 
tion?. Superimposing the excluded region from 
the day-night spectra in Fig. 12 with the allowed 
region from the global rates in Fig. 9 we see that 
almost the full allowed region at 95 %CL lays be- 
low the 95% CL excluded region from the day- 
night spectra. However this does not mean that 
the SMA solution is ruled out at that CL. The re- 
sult from the correct statistically combined analy- 
sis is what we show in Fig. 13 and more important 
in Table 1. From this results we learn that the 
SMA can describe the full data set with a prob- 
ability of 34%, but it is now shifted to smaller 
mixing angles to account for the flatter spectrum. 

(d) Similar statement as (c) holds for the SMA 
solution for sterile neutrinos. 

Thus our conclusion is that from the statisti- 
cal point of view all solutions are acceptable since 
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they all provide a reasonable GOF to the full data 
set. LMA and LOW-QVO solutions for oscilla- 
tions into active neutrino seem slightly favoured 
over SMA solutions for oscillations into active or 
sterile neutrinos but these last two are not ruled 
out. 

Let's comment now how these conclusions may 
depend on the specific features of the SSM model 
and/or the analyzed data. In order to do so we 
are going to discuss two departures from the pre- 
vious global analysis: 

(a) effect of the uncertainty in the predicted 
SSM boron flux, and 

(b) effect of removing from the analyzed data 
one of the measured rates. 

Concerning (a) in Fig. 14 we show the results 
of the analysis for active and sterile oscillations 
when we allow the departure from SSM normal- 
ization for the 8 B flux. In doing so we have 
treated the boron flux as a free parameter that 
must be fixed by the experiments, in particular by 
SK. The quality of the different minima is shown 
in the last line of Table 1. Comparing with the 
results for fixed SSM 8 B flux we see that allow- 
ing free S B normalization leads to both an im- 
provement of the quality of the SMA and LMA 
solutions but the effect is more important for the 
SMA. The small improvement of the LMA solu- 
tion is due to the fact that allowing for a free 
S B flux leads to a better simultaneous agreement 
with CI and SK measured total rates. The most 
important effect on the SMA solution arises from 
the fact that in order to account for a flatter spec- 
trum the SMA solution has shifted to smaller an- 
gles, what also implies a larger S B flux contribu- 
tion to both CI and SK. This can be compensated 
by the allowance of a lower predicted 8 B flux. 
Notice also that the LOW solution and QVO so- 
lutions which have the best fit very close to the 
SSM 8 B normalization, the quality of the fit has 
decreased simply due to the effect of the addi- 
tional free parameter. 

For (b) the situation is more involved. In 
Fig. 15 we show the analysis for active and sterile 
oscillations where we have removed the number 
of events measured by each experiment, SK, CI 
and Ga respectively, from the global analysis but 
always keeping the SK day-night spectra. The 



quality of the local minima is showed in Table 2. 
From the figure and the table we can draw two 
main conclusions: 

- The results for oscillations into active and ster- 
ile neutrinos are only significantly different for the 
LMA and LOW-QVO solutions when both SK 
and CI rates are included in the analysis. This is 
consequence of the larger survival probability for 
S B neutrinos needed to accommodate the mea- 
sured rate at SK for oscillations into sterile neu- 
trino (due to the absence of the NC contribu- 
tion) which for large mixing angles lead to a too 
large contribution to Chlorine rate. On the con- 
trary, if the chlorine experiment is removed from 
the global analysis both oscillations into active or 
sterile neutrinos provide equivalently good fits to 
the data and in particular a large region in both 
sides of maximal mixing is allowed [17] . 

- The SMA region is disfavoured in the global 
analysis due to the spectrum of SK and this re- 
sult holds independently of the removal of one 
the measured rates. It is clear by comparing the 
active SMA and LMA columns of Table 2 that 
the relative quality of SMA versus LMA worsens 
when SK or CI total rates are removed from the 
analysis. However once both are included (third 
row), the difference between SMA and LMA de- 
scriptions decreases because, as explained before, 
the worse SMA fit to the SK spectrum is compen- 
sated by the better fit to these two rates simul- 
taneously As a matter of fact, in this case the 
best solution lies in the LOW region where both 
rates can be better fitted than for LMA and the 
spectrum data is better described than for SMA. 

Finally let's stress that when the allowed pa- 
rameter space for solutions of a physics problem 
consists of a set of isolated regions, the quality of 
a given solution cannot be inferred from the size 
of the corresponding region. This fact is nicely il- 
lustrated when comparing the SMA solutions for 
the active and sterile cases shown in Fig. 13. No- 
tice that the SMA region for oscillations into ster- 
ile neutrinos is larger than the corresponding one 
for active neutrinos. This is because it is defined 
with respect to the global minimum for the ster- 
ile oscillation hypothesis (which lies in the SMA 
region) while the SMA region for the active case 
is defined with respect to the global minimum 
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Figure 12. 99 (shadow), 95 and 90 excluded regions from the analysis of the SK measured spectrum at 
night and at day. The inverted triangle represents the minimum used to define the regions. 
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Figure 13. 90, 95 and 99 % CL allowed regions from the global analysis of solar neutrino data including 
the total measured rates and the SK measured spectrum at day and night. The global minimum is 
marked with a star while the local minima are denoted with a dot. 
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Table 1 

Best fit points and GOF for the allowed solutions for different combinations of observables. 



Observable 




Active 


Sterile 


Rates 




SMA 


LMA 


LOW 


VAC-QVO 


SMA 


Am'/eV z 
tan 2 9 


5.5 x 10" b 
0.0015 


1.9 x lO" 5 
0.3 


9.2 x 10" 8 
0.65 


9.7 x 10- 11 
0.51 (1.94) 


4.1 x 10" b 
0.0015 


Prob (%) 


50% 


8% 


0.5% 


2% 


19% 


Rates 
+Zcnith 


Am 2 /eV 2 
tan 2 6» 


5.6 x 10" b 
0.0012 


5.2 x 10" b 
0.35 


7.9 x 10" 8 
0.69 


9.7 x lO" 11 
0.51 (1.94) 


4.1 x 10" b 
0.0015 


Rates 
+Spectrum 


Am 2 /eV 2 
tan 2 9 


5.0 x 10" b 
0.00075 


2.1 x 10" 5 
0.32 


9. x 10" 8 
0.67 


8.4 x i0~ lu 
1.7 (QVO) 


3.9 x 10-° 
0.00069 


Rates 
+Spec£i+Spec]v 


Am 2 /eV 2 
tan 2 9 


5.0 x 10 _b 
0.00058 


3.2 x 10" 5 
0.33 


1. x 10"' 
0.67 


8.6 x lO" 1 " 
1.5 (QVO) 


3.9 x 10" b 
0.0006 


Prob (%) 


34% 


59% 


40% 


29% 


30% 


Rates 

+Spec£i+Specjv 

free S B 


Am 2 /eV 2 
tan 2 9 


4.9 x 10" b 
0.00046 


3.0 x 10" 5 
0.26 


9.7 x 10" 8 
0.67 


8.4 x nr iU 

1.7 (QVO) 


3.9 x 10" b 
0.00053 


Prob (%) 


51 % 


65% 


36% 


27% 


42% 
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Figure 14. 90, 95 and 99 % CL allowed regions from the global analysis of solar neutrino data taking the 
boron flux as a free parameter. The global minimum is marked with a star while the local minima are 
denoted with a dot. 
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Table 2 

Xmin an d GOF in the regions of SMA, LMA and LOW for the global analysis when one of the fluxes is 

removed from the data (35 dof). 



Observable 




Active 


Sterile 


Rates (Ga + CI) 
+SpecD +SpecAr 




SMA 


LMA 


LOW 


SMA 


LMA 


LOW 


Xmin 

Prob (%) 


38.9 
30% 


29.8 

72% 


34.3 

50% 


39.2 

29% 


29.6 

73% 


34.3 

50% 


Rates (Ga + SK) 
+SpecD +Specjv 


/Cmin 
Prob (%) 


38.1 

33% 


30.5 
69% 


30.4 
69% 


37.2 

37% 


31.4 

64% 


29.4 

73% 


Rates (CI + SK) 
+SpecD +Spec./v 


\rnin 

Prob (%) 


38.0 

33% 


33.0 

56% 


31.3 

65% 


39.7 

27% 


45.7 
12% 


37.8 

34% 



for that scenario which lies in the LMA. However 
when looking at the GOF of these two solutions 
we find that the SMA solution for active case pro- 
vides a better description of the data than the 
sterile one despite its region is smaller. In order 
to put the two scenarios, active and sterile in the 
same footing we must embed them into a common 
framework as we describe in the next section. 



4. Four Neutrino Oscillations 

In this section we present a brief update of the 
analysis performed in Ref. [22]. We refer to this 
publication for further details on the as well as 
for the relevant references. Here we simply sum- 
marize the main ingredients. 

Together with the results from the solar neu- 
trino experiments we have two more evidences 
pointing out towards the existence of neutrino 
masses and mixing: the atmospheric neutrino 
data and the LSND results. All these experi- 
mental results can be accommodated in a single 
neutrino oscillation framework only if there are 
at least three different scales of neutrino mass- 
squared differences. The simplest case of three in- 
dependent mass-squared differences requires the 
existence of a light sterile neutrino, i.e. one whose 
interaction with standard model particles is much 
weaker than the SM weak interaction, so it does 
not affect the invisible Z decay width, precisely 
measured at LEP. 

In four-neutrino schemes the flavor neutrino 
fields v a L (a — e, s, /i, t) are related to the fields 



VkL of neutrinos with masses mk by the relation 

4 



V a L 



E 
fe=i 



UakVkL (a = e,s,fj,,T), 



(13) 



where U is a 4x4 unitary mixing matrix, which 
contains, in general 6 mixing angles (we neglect 
here the CP phases). 

Existing bounds from negative searches for 
neutrino oscillations performed at collider as well 
as reactor experiments impose severe constrains 
on the possible mass hierarchies as well as mix- 
ing structures for the four neutrino scenario. In 
particular they imply: 

(a) Only two four-neutrino schemes can ac- 
commodate the results of all neutrino oscilla- 
tion experiments. In both these mass spectra 



A 



B 



solar 



atm 



t 
LSND 



t 
LSND 

I 



atm 



solar 



there arc two pairs of close masses separated 
by a gap of about 1 eV which gives the mass- 
squared difference Am 2 SBL = Am^ responsible 
for the short-baseline (SBL) oscillations observed 
in the LSND experiment (we use the common 
notation AmL = m\ — m3). We have ordered 
the masses in such a way that in both schemes 

Am L« = Am 
tions and Am 



2 

21 
2 
atm 



produces solar neutrino oscilla- 
is responsible for at- 



'•43 



mospheric neutrino oscillations. 
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(b) In the study of solar neutrino oscillations 
only four mixing angles are relevant and the U 
matrix can be written as 



U — U34 U24 U23 U\2 ■ 

where #12, $23, ^24, 
and we will define a 



(14) 

$ 34 are four mixing angles 
l3 = cos ■dij and sy = sin$jj. 

Since solar neutrino oscillations are generated 
by the mass-square difference between V2 and v\ , 
it is clear from Eq. (14) that the survival of so- 
lar z/ e 's mainly depends on the mixing angle $12, 
whereas the mixing angles #23 and $24 determine 
the relative amount of transitions into sterile v s 
or active v^ and v T . Let us remind the reader 
that v^ and v T cannot be distinguished in solar 
neutrino experiments, because their matter po- 
tential and their interaction in the detectors are 
equal, due only to NC weak interactions. The ac- 
tive/sterile ratio and solar neutrino oscillations in 
general do not depend on the mixing angle $34, 
that contribute only to the different mixings of 
Vp and v r , and depends on the mixing angles $23 
$24 only through the combination cos $23 cos $24- 

Therefore, the oscillations of solar neutrinos de- 
pend only on $12 and the product cos $23 cos $24- 
If cos $23 cos $24 7^ 1, solar VeS can transform in 
the linear combination v a of active v^ and v T . We 
distinguish the following limiting cases: 

• cos $23 cos $24 = corresponding to the limit 
of pure two-generation v e — > v a transitions. 

• cos $23 cos $24 = 1 for which we have the limit 
of pure two-generation v e — ► v s transitions. 

In the general case of simultaneous v e — ► v s 
and v e — > v a oscillations the corresponding prob- 
abilities are given by 



P., 
P v 



-23^24 



(1 - P V . 



= (1 - 4 3 C 2 24 ) (1 - P„.^„.) 



(15) 
(16) 



where P Ve ^ Ve takes the standard two-neutrino 
oscillation form (3) for Am^ and 0\2 but com- 
puted for the modified matter potential 



A = A C c + cl 3 cl A A NC ■ 



(17) 



Thus the analysis of the solar neutrino data in 
the four-neutrino mixing schemes it is equiva- 
lent to the two-neutrino analysis but taking into 
account that the parameter space is now three- 
dimensional (Am|i, tan 2 $12, cos 2 $23 cos 2 $24)- 



We first present the results of the allowed re- 
gions in the three-parameter space for the global 
combination of observables. In building these re- 
gions, for a given set of observables, we com- 
pute for any point in the parameter space of 
four-neutrino oscillations the expected values of 
the observables and with those and the corre- 
sponding uncertainties we construct the function 
X 2 (Ato 2 17 $12, C23C24). We find its minimum in 
the full three-dimensional space. The allowed re- 
gions for a given CL are then defined as the set 
of points satisfying the condition 
X 2 (A™ 2 1 ,$ 1 2,c 2 3 c 2 4 ) - x 2 min < Ax 2 (CL, 3 dof) 
where, for instance, Ax 2 (CL, 3 dof)=6.25, 7.83, 
and 11.36 for CL=90, 95, and 99 % respectively. 
In Figs. 16 we plot the sections of such volume in 
the plane (Am|i, tan 2 ($i2)) for different values of 
c 23 c 24- The global minimum used in the construc- 
tion of the regions lies in the LMA region and for 
pure active oscillations value of C2 3 c 2 4 = 0. 

As seen in Fig. 16 the SMA region is always a 
valid solution for any value of C2 3 c 2 4 . This is ex- 
pected as in the two-neutrino oscillation picture 
this solution holds both for pure active-active 
and pure active-sterile oscillations. Notice, how- 
ever, that the statistical analysis is different: in 
the two-neutrino picture the pure active-active 
and active-sterile cases are analyzed separately, 
whereas in the four-neutrino picture they are 
taken into account simultaneously in a consistent 
scheme. We see that in this "unified" framework, 
since the GOF of the SMA solution for pure ster- 
ile oscillations is worse than for SMA pure ac- 
tive oscillations (as discussed in the previous sec- 
tion) , the corresponding allowed region is smaller 
as they are now defined with respect to a common 
minimum. 

On the other hand, the LMA, LOW and QVO 
solutions disappear for increasing values of the 
mixing c| 3 C2 4 . In Fig. 17 we show the the val- 
ues of A% 2 (with respect to the global minimum) 
for the different solutions as a function of c 2 3 C2 4 . 
From this figure we can read the limiting values 
of c 2 3 C2 4 for which a given solution is allowed in 
a given CL. We list some of this limits in Table 3. 
We find that at 95% CL the LMA is allowed for 
maximal active-sterile mixing c\ 3 c\^ = 0.5 while 
at 99% CL all solutions are possible for this max- 



1G 




tan 2 (i? 12 ) 



Figure 16. Results of the global analysis for the allowed regions in Ara^ and tan 2 $12 for the four- 
neutrino oscillations. The different panels represent the allowed regions at 99% (darker) and 90% CL 
(lighter) obtained as sections for fixed values of the mixing angles c\^c^ A of the three-dimensional volume. 
The best-fit point in the three parameter space is plotted as a star. 
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Figure 17. A% 2 as a function of the mixing pa- 
rameter c| 3 C24 for the different solutions SMA 
(full line), LMA (dashed), LOW(dotted) and 
QVO (dot-dashed) from the global analysis of the 
rates and the day-night spectrum data. The dot- 
ted horizontal lines correspond to the 90%, 95 %, 
99% CL. 



Table 3 

Maximum allowed value of COS23 cos 24 a * 90%, 
95%, and 99 % CL for the different solutions to 
the solar neutrino problem. 

CL SMA LMA LOW QVO 



90 


0.9 


0.44 


0.3 


forbidden 


95 


all 


0.53 


0.44 


0.28 


99 


all 


0.72 


0.77 


0.88 



imal mixing case. 

Other authors have recently performed fits of 
the atmospheric neutrino data in the framework 
of the two four-neutrino schemes A and B [23]. 
We expect that in the future a combined fit of 
solar and atmospheric neutrino data will allow to 
constraint further the mixing of four neutrinos. 

5. Summary and Conclusions 

We have studied the status of the solutions of 
the solar neutrino problem in terms of oscillations 
of v e into active or sterile neutrinos in the frame- 
work of two-neutrino and four-neutrino mixing, 
after a global analysis of the full data set cor- 
responding to the latest SK data on the total 



event rate, their day-night dependence and the 
recoil electron energy spectrum, together with the 
data from Chlorine and Gallium experiments pre- 
sented at the ^-2000 conference. In particular 
we have shown the possible solutions in the full 
parameter space for oscillations including both 
MSW and vacuum, as well as QVO and matter 
effects for mixing angles in the second octant 

Our main conclusions from the two-neutrino 
oscillation analysis are the following: 

(a) From the statistical point of view, all solu- 
tions for oscillations into active neutrinos: 
LMA, LOW, SMA, QVO solutions are ac- 
ceptable since they all provide a reasonable 
GOF to the full data set. 

(b) The same holds for the SMA solution, which 
is the best possible solution for oscillations 
into sterile neutrinos. 

(c) LMA and LOW-QVO solutions for oscil- 
lations into active neutrinos seem slightly 
favoured over SMA solutions for oscillations 
into active or sterile neutrinos due to the 
flatter spectrum of SK. 

(d) Oscillations into sterile neutrino provide a 
slightly worse fit. 

(e) When allowing for a free 8 B normalization 
weakens the conclusion (c). 

(f ) Removing the total rate measured at Chlo- 
rine or SK from the analysis weakens con- 
clusion (d) and straightens conclusion (c). 

For the analysis of four-neutrino oscillations, 
which allows for oscillations into a state which is 
a combination of active and sterile neutrino, we 
find that the SMA region is always a valid solu- 
tion while the LMA, LOW and QVO solutions 
disappear for increasing values of the additional 
mixing C23C24. However we find that at 95% CL 
the LMA is still allowed for maximal active-sterile 
mixing C23C24 = 0.5 while at 99% CL all solutions 
are possible for this maximal mixing case. 
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